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Q\ • In lyj, a new method was developed to calculate energy levels in 

QED nonrelativistic bound states, up to order ma^ (or a^ with re- 
Qh! spect to the Bohr levels). Whether or not this method could be used 

o ' beyond this order was left as an open question. In this paper, we 

1^ ■ answer this question with the help of a nonrelativistic effective field 

theory, NRQED. This theory permits to deal separately with bound 

state physics (taking place at momenta of order ma or below) from rel- 

S^ ■ ativistic physics (where the momentum is of order the electron mass) 

in a rigorous and systematic manner. We find that, aside from in- 
frared terms which had to be discarded without justification in ||l|, 
the 0{ma^) and 0{ma^) calculations give the same results in both 
methods. It is shown, however, that this is due to the special physical 
origin of these contributions and that the method of Q would not 
give the right answer at order ma^ or higher. The separation of scales 
provided by NRQED is essential to the derivation of this result. 

1 Introduction 

In III] , fine and hyperfine splittings calculations of nonrelativistic QED bound 
states are carried out up to order ma^, or a^ with respect to the lowest order 



(Bohr) energy levels. The shifts in energy are obtained, essentially, by sand- 
wiching Feynman scattering amplitudes between Schrodinger wavef unctions. 
These amplitudes are evaluated with the external particles on mass-shell and 
their four- momenta taken in the nonrelativistic limit. The authors left open 
the question of whether this approach could be pushed beyond this order. 

In this paper, we answer this question using a nonrelativistic effective 
field theory, nonrelativistic quantum electrodynamics (NRQED)|0]. The use 
of such an effective theory is possible because the typical QED bound state 
momentum is of order fia (where fi is the reduced mass), which is much 
smaller than the masses of the constituents. In an effective field theory, an 
infinite number of nonrenormalizable interactions are present, but they are 
suppressed by factors of 1/A where A is an appropriate cutoff characterizing 
a scale of energy at which the theory would break down. Because of the 
nonrelativistic nature of the bound states we are considering, it is natural to 
take this cutoff as being the reduced mass of the system. It is important to 
emphasize that the physics taking place at energies beyond the cutoff A does 
affect the low energy physics and thus must not be thrown away. However, 
because this physics takes place at such a high energy (in comparison to the 
domain of energy in which the effective theory is used), it can be incorpo- 
rated back in the low energy theory by renormalizing its coefficients. Let us 
emphasize that this procedure is performed in a rigorous and systematic way 
and that the low energy theory can thus be made to reproduce QED to an 
arbitrary precision (as long as the system under study is characterized by 
momenta p < A). 

We find that the approach of ^ would be inconsistent at the C(ma^) and 
would not lead to the correct shift in energy. Moreover, using NRQED, we 
will be able to understand why this approach gives the correct order (ma^) 
and ma^ contributions^. 

This paper is organized as follows. In Section 1, we explain the difficulties 
associated with QED nonrelativistic bound states. In the following two sec- 
tions, we show how to build a nonrelativistic effective field theory and why 
it is particularly useful for the study of bound states. We then evaluate the 
hyperfine splitting in Positronium up to order ma^ and show what distin- 

^Except for the infrared divergent terms which, in the formalism of ||l|, must be set 
aside without any further justification. In the context of NRQED, however, these terms 
are systematically removed. 



guishes the 0{a'^) and 0{a^) terms from higher order corrections. This will 
clearly show why the approach of [|^ is valid at that order (except for the 
infrared cutoff dependence), and also why one can generalize the order ma^ 
result to arbitrary values of n, the principal quantum number. We conclude 
by explaining how to proceed to obtain the ma^ and higher corrections, and 
why the approach of Ref . [|l| would fail at this stage. 

2 Nonrelativistic bound states 

In a nonrelativistic bound state, the traditional expansion in the number 
of loops familiar from quantum electrodynamics (QED) fails. The physical 
reason is the presence of energy scales absent from scattering theory. Indeed, 
the finite size of the atom provides an energy scale of order the Bohr radius 
Ri 1/ma which, by the uncertainty principle, is equal to the inverse of the 
typical three-momentum of the constituent particle (~ ma). Because of 
this new energy scale, there is a region of the momentum integration in 
which the addition of loops will not result in additional factors of a. Also, 
because this energy scale is much smaller than the particles masses, the 
system is predominantly nonrelativistic, a simplification not taken advantage 
of in traditional approaches. Moreover, a third energy scale, again vastly 
different from the previous two, is set by the particles kinetic energies and 
further complicates bound state calculations. It is at the origin of the so- 
called retardation corrections which lead, for example, to the Lamb shift. 

To illustrate more specifically the problems arising from the presence of 
these energy scales, consider a traditional Bethe-Salpeterp] calculation of the 
energy levels of positronium. The potential in the Bethe-Salpeter equation is 
given by a sum of two-particle irreducible Feynman amplitudes. One gener- 
ally solves the problem for some approximate potential and then incorporates 
corrections using time-independent perturbation theory. Unfortunately, per- 
turbation theory for a bound state is far more complicated than perturbation 
theory for, say, the electron (^-factor. In the latter case a diagram with three 
photons contributes only to order a^. In positronium a kernel involving the 
exchange of three photons contributes not only to order a^, but also to all 
higher orders as wellQ: 

^We assume here that the Coulomb com.ponent of each photon exchange has been pro- 
jected out. See Section 4 for more details on the special status of the Coulomb interaction 



(Vi) = a^m(ao + aia + 020^ + ...)• (1) 

Therefore, in the bound state calculation there is no simple correlation be- 
tween the importance of an amplitude an the number of photons in it. 
Such behavior is at the root of the complexities in high-precision analyses of 
positronium or other QED bound states, and it is a direct consequence of 
the multiple scales in the problem. Any expectation value like Eq.(l) will be 
some complicated function of the ratios of the three scales in the atom: 

(1/i) = a^mF{{p)/m,{K)/m). (2) 

Since {p)/m ~ a and {K)/m ~ a^, a Taylor expansion of F in powers of 

these ratios generates an infinite series of contributions just as in Eq.(l). A 
similar series does not occur in the gf-factor calculation because there is but 
one scale in that problem, the mass of the electron. 

Let us again emphasize that traditional methods for analyzing these 
bound states fail to take advantage of the nonrelativistic character of these 
systems; and atoms like positronium are very nonrelativistic: the probability 
for finding relative momenta of 0{m) or larger is roughly a^ c:^ 10^^^! 

Let us also point out that the difficulties associated with the many dif- 
ferent energy scales in bound states are not limited to QED. Consider the 
study of nonrelativistic QCD bound states, the T for example. Even if the 
means of extracting the properties of such systems, namely simulating the 
system on a lattice, are very different than the use of the Bethe-Salpeter 
equation, the nonrelativistic nature of the bound state still greatly compli- 
cates calculations. Indeed, the space-time grid used in such a simulation 
must accommodate wavelenghts covering all of the scales in the meson, rang- 
ing from {mv'^)~^ down to m~^. Given that f^ ~ 0.1 in T, one might easily 
need a lattice as large as 100 sites on side to do a good job. Such a lattice 
would be three times larger than the largest wavelenght, with a grid spacing 
three times smaller than the smallest wavelenght, thereby limiting the er- 
rors caused by the grid. This is a fantastically large lattice by contemporary 
standards and is quite impractical. The strong force equivalent of NRQED, 
NRQCD0], greatly reduces the size of the lattices required to study heavy 
quarks systems since it essentially eliminates the mass of the meson as a 
relevant energy scale in the calculation (we will soon see how this occurs in 
NRQED). 

in nonrelativistic bound states. 



3 NRQED 

To take advantage of the simplifications associated with nonrelativistic sys- 
tems, one can approximate QED in the hmit p/m ^ 1 (which is equivalent 
to the limit v/c ^1). One way to do this consists of simply expanding the 
QED Lagrangian in powers oi p/m to obtain NRQED. The usefulness of this 
expansion is enhanced by performing a Foldy-Wouthuysen transformation, 
which decouples the electron and positron degrees of freedom. It is also, how- 
ever, possible to define NQRED without prior knowledge of QED. Indeed, 
one can build up the Lagrangian of NRQED by imposing restrictions coming 
from the symmetry obeyed by the theory, such as gauge invariance, chiral 
symmetry in the limit mg -^ 0, and Lorentz invariance for the photon kinetic 
term. Lorentz invariance for the rest of the Lagrangian and renormalizability 
are not necessary. 

It is of course possible to write down an infinite number of terms respect- 
ing these symmetries, but this is not a problem since operators of dimension 
greater than 4 will be multiplied by coefficients containing inverse powers of 
m, the only available energy scale in the problem (remember that m plays 
the role of the ultraviolet cutoff A in NRQED). Since we are considering the 
limit p/m <^ 1, these operators will yield contributions suppressed by that 
many powers of p/m and therefore, for a given accuracy, only a few terms 
need to be kept in an actual calculation. The first few terms one obtains are 
then 

-Cnrqed = -l/A{F^''f + ^lj^i^idt-eAo + B^/2m + Tiy8m^ 
+ci e /2m cr ■ B + C2e /8m^ V • E 
+C3 ie/8m^cr- (D x E - E x T>) + ...\^ 

+die^ /m^ (^V)' + ^e' M' (^V^)^ + ... 

-|- positron and positron-electron terms. (3) 

where the c's and d's are coefficients of order one. Here, ip and x ^-^e the 
(two-component) electron and positron spinors and D is the gauge-invariant 
derivative. An example of a positron-electron term (and the one we will focus 
on later) is d^e^/m"^ (ip'^crx) ' (x^'^'V')) which represents, to lowest order, the 
process e~^e~ — > e~^e~ in the s channel. This term, and the other two having 
coefficients starting with rf's are examples of "contact" interactions. 



Notice that we have recovered the interactions famihar from nonrelativis- 
tic quantum mechanics such as the first order relativistic correction to the 
kinetic energy (D^/8m^) and the Darwin term (c2 e/8m^ V • E). 

Of course, one cannot perform a nonrelativistic expansion of the photon 
kinetic term. It is however still possible to simplify the nonrelativistic calcu- 
lation by choosing the most useful gauge in that context, the Coulomb gauge 
(the reason for its usefulness will become clear later on). In that gauge, the 
Coulomb propagator is given by i/k^ and the transverse propagator is given 
by 

The non recoil limit is obtained by setting /c^ = — k^ in the transverse pho- 
ton propagator, which amounts to neglecting the energy transferred by the 
photon. 

We now have to fix the coefficients appearing in £nrqed- To do so, we may 
simply compare scattering amplitudes in both QED and NRQED, at a given 
kinematic point (which we choose to be at threshold, i.e. with the external 
particles at rest). This renormalization of the low energy coefficients (often 
referred to as the "matching procedure" in the literature on effective field 
theories) must be carried out to as many loops as required for the accuracy 
desired in the final calculation; at that stage, one can associate a factor of e for 
each vertex to evaluate the contribution to the bound state calculation (see 
section [4.1] for more details). Only a finite number of NRQED interactions, 
for that given number of loops, must be kept. That number is determined 
using NRQED counting rules, which permit to determine the order at which 
a given NRQED diagram will contribute, a/ter renormalization (or matching) 
has been carried out. We will give these rules in section [4.1]. 

4 Bound states diagrams in NRQED 

Now that NRQED is completely defined, we can apply it to bound state cal- 
culations, where its usefulness is most apparent. To be specific, consider the 
energy levels of positronium. In a bound state, the expansion in Feynman 
diagrams of NRQED reduces to conventional Rayleigh-Schrodinger pertur- 
bation theory. The "external" wavefunctions of course don't correspond to 
free wavepackets but to bound states. In the traditional Bethe-Salpeter for- 
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inalisin, the form of these wavefunctions depend on which interactions one 
desires to include exactly, and which interactions one wants to treat per- 
turbatively. Depending on this choice, the external wavefunctions may be 
as simple as the Schrodinger wavefunctions or as complicated as the Dirac 
wavefunctions. 

We are now in measure to appreciate the usefulness of the Coulomb gauge. 
This gauge permits us to keep the simplest interaction for the zeroth order 
kernel, namely the Coulomb interaction, in which case the Bethe-Salpeter 
equation reduces to the Schrodinger equation and the wavefunctions are the 
ones found in any textbook on quantum mechanics. In this calculation we will 
work only with the IS" wavefunction in positronium (reduced mass /i = m/2) 
which is given, in momentum space, by 



(p2 + ^2 



*(P) = T-TT^., (5) 



where 7 is the typical bound state momentum, equal to fiv = ma/2 in 
positronium. The ground state energy is — 7^/2/i = —ma'^/A. Physically, an 
infinite number of Coulomb interactions is incorporated in the wavefunction. 
The reason why at least an infinite number of Coulomb kernels must be 
taken into account in the wavefunction can be easily understood by con- 
sidering fig. [la], where a generic bound state diagram is shown. Ki and 
K2 represent some unspecified kernels. We can schematically represent this 
diagram by 

^'^F(k;k) (6) 



(27r)3' 

where F(pi;p2) is a function including the bound state wavefunctions, the 
interactions present in the kernels Ki and K2, and the propagator of the 
electron-positron pair between these two kernels. The arguments pi and 
P2 represent the three-momenta flowing out of the first kernel and into the 
second one. In figure [la], Pi = P2 = k. As explained above, working with 
the effective field theory means that the electron mass no longer appears as a 
dynamical energy scale in the NRQED diagrams; this means that m can only 
appear as an overall factor of the form l/m" in figure [la]Q. The function F is 
therefore independent of the electron mass; the only energy scale it depends 

•^In the language of eflFective field theories (EFT's), NRQED is a decoupHng EFT with 
cutoff scale A = m. The heavy particle being decoupled is the electron's antiparticle so 
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on is the typical bound state momentum 7, and the ultraviolet cutoff (when 
it is divergent )0. 

Let us now consider adding a Coulomb interaction between the two ker- 
nels, as in fig. [lb]. In conventional scattering theory, one would conclude 
that this diagram is of higher order than diagram [la] because of the fac- 
tor a coming from the vertices. However, this argument does not hold in a 
bound state, because the typical momentum flowing through the Coulomb 
interaction is of order 7. Indeed, consider how the integral Eq. [^ is modified 
when one adds the Coulomb interaction. It must now be written 

d^k(fl ie^ i 

l2nf |k-l|2 (_y-k2)/m ^ ' '' 

/■ Skd^l 1 

= '""^" i (2vr)e(y + k2)(|k-l|2) -^(^' ^) (^) 

where i/\k — Ip is the Coulomb propagator and i?7i/(— 7^ — k^) is the non- 
relativistic propagator for the e~e^ pair. 

We see that the mass factors out (more about this later!), leaving 7 as the 
only energy scale in the integral. By dimensional analysis, the final result will 
be of order ~ ma (the overall factor) divided by /7 (from the integral) ~ 1. 
This shows that in a bound state, adding a Coulomb interaction does not, 
in fact, lead to a diagram of higher order. This is why one must sum up all 
Coulomb interactions in the wavefunctions right from the start. This implies 
the relation expressed in fig. [2] which is nothing else than a diagrammatic 
representation of the Schrodinger equation. 

4.1 Counting rules 

The argument we used to show that adding a Coulomb interaction does not 
increase the order of a bound state diagram is an example of a NRQED 

that Lorentz invariance is broken (as expected in a nonrelativistic theory). Of course, 
there is a positron in NRQED, but it is here an independent particle; it never appears 
as the negative energy pole of the electron propagator (in the language of time-ordered 
perturbation theory, there are no Z graphs). Consequently, one can apply the decoupling 
theorem g to NRQED, with the consequence that the heavy scale m will only appear as 
an overall factor l/m" in the interactions. One also expects m to enter in the running of 
the low energy coefficients through logarithms, but this dependence on m will only come 
from the QED scattering diagrams, through the matching (see section [5.3]. 
"^We will deal with the infrared divergences separately. 



counting rule. These rules are extremely important since they permit one to 
estimate the order of contribution of a diagram without calculating it. Let us 
first consider non-recoil diagram^. In that approximation, the counting rules 
are based on the observation that the mass m of the electron always factors 
out of the integrals so that it no longer represents a relevant energy scale. 
This is clear since m factors out of the nonrelativistic fermion propagators, 
as we saw in Eq.[0, and it enters the rules of the NRQED vertices only as 
an overall factor. We therefore see that the use of NRQED eliminates the 
relativistic energy scale of the system, which, in the light of the discussion 
in the introduction, greatly simplifies the analysis. 

Starting from this observation, we can recast the discussion of the Coulomb 
interaction in the following manner: adding a Coulomb interaction between 
two kernels as in fig.[l] leads to an additional factor of a coming from the 
vertices and a factor of m coming from the additional fermion propagator. 
The corresponding overall factor of ma must be canceled by a factor having 
the dimensions of energy in order to keep the dimensions straight, but since 
the only energy scale left in the problem is 7, we finally obtain a correction 
of order ma/7 — 1- It now becomes obvious that the Coulomb interaction 
is the only kernel having the property of not increasing the order of a bound 
state diagram since all the other interactions contain additional factors of 
1/m. Each of these 1/m factors will have to be canceled by a corresponding 
factor of 7, leading therefore to a result of higher order in a. 

With these rules, it becomes extremely simple to evaluate the order of the 
contribution of a given bound state diagram. All we need to know is that the 
external wavefunctions contribute a factor m^a^. To obtain the contribution 
of a given kernel, one has simply to count the number of explicit factors of a 
(coming from the vertices) and inverse factors of m (coming from the vertices 
and the e~e~^ propagators) this diagram contains. After canceling the 1/m's 
with factors of 7, the number of a's left over give the order of the contribu- 
tion of the diagram. Consider, for example, inserting a transverse photon (in 
the non-recoil limit) between two wavefunctions as in fig. [3]. The transverse 

^In the case of a bound state having equal mass constituents, as positronium, there 
is no distinction possible between "recoil effects" and the so-called "retardation effects" . 
This is not, for example, the case in muonium or hydrogen where (sometimes confusing) 
distinctions are made between non-recoil and non-retardation approximations. What we 
mean here by the non-recoil approximation is the limit in which the transverse photon 
does not carry any energy. 



photon vertex has one power of e and one power of 1/m. Taking the two 
vertices into account, we can conclude that this diagram will contribute to 
order m?a^ x a/m? ~ ma^ (it is not necessary to multiply by any factor of 
7 here). The same is true of the annihilation diagram and of the spin-spin 
interaction as they also contain a factor of 1/m? relative to the Coulomb 
interaction. Another important observation is that these diagrams will con- 
tribute to only one ordeij^ in a, in contradistinction with Feynman diagrams 
in conventional Bethe-Salpeter analysis. It is important to understand that 
these rules give the order of contribution of a diagram after the renormal- 
ization of the effective theory has been performed. Indeed, many NRQED 
bound state diagrams are actually badly divergent in the ultraviolet. The 
rules give us the order of the finite contributions left over after the NRQED 
bare coefficients have been renormalized (or, in other words, after matching 
with QED has been performed). 

The situation is only slightly more complicated when one takes into ac- 
count recoil effects. We will not dwell on this issue here, but let us just 
mention that they make the kinetic energies K enter as a mass scale and 
that they are at the origin of the appearance of log's of a which are charac- 
teristic of bound state^. 

Now that we have a way to estimate the order (in powers of a) of the 
NRQED diagrams contributions, we can go back to the matching of the ef- 
fective Lagrangian (Eq. [Q) coefficients with QED. As mentioned earlier, this 
matching is done by evaluating scattering amplitudes evaluated in both the- 
ories. There are two expansions involved in this matching: an expansion in 
number of loops in both QED and NRQED, and an expansion in p/m in 
NRQED. The number of loops n required is found by simply setting ma^^"^ 
equal to the precision desired in the final result (a factor of m?a^ comes from 
the wavefunctions and a factor of l/rn^ comes from the nonrelativistic nor- 
malization of the Feynman amplitudes). For a given number of loops, there 
is an infinite number of nonrelativistic interactions in the effective theory 



^We are still limiting ourselves to the non-recoil limit. 

^To be more precise, the logarithms involve the bound state typical energy scale 7 ~ ma 
and the heavy scale A = m, leading to a result of the form \n{malm) = In a. This log 
of m is the one we mentioned earlier when talking about the running of the low energy 
coefficients. Notice that the infrared cutoff A always gets canceled while performing the 
matching because, by construction, the infrared behavior of the effective theory is the 
same as the one of the full theory. 
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but, using the counting rules described above, only a finite number must be 
kept for a given accuracy in the final result. Even if the rules were derived 
for bound state diagrams, they can be directly applied to the NRQED scat- 
tering diagrams i.e., one imagines sandwiching them between wavefunctions 
and one performs the same power counting. Once again there are strong 
ultraviolet divergences but let us repeat that the counting rules give us the 
order at which the diagrams will contribute once the matching of the NRQED 
coefficients will have been performed (in other words, once the counterterms 
will have been included). 

There is only one exception to the counting rules, and it occurs when 
a NRQED scattering diagram contains a power-law infrared divergence. In 
that case, the counting rules must be modified to take into account the pres- 
ence of a new scale, the infrared cutoff. Take, for example, a scattering 
diagram which has an overall factor of m'^a^ times an integral J d^kf{]<.,'j). 
The counting rules would suggest that this diagram will contribute by a term 
of order m'^a^/'j = ma^ to the bound state calculation but, if it is linearly 
infrared divergent, it can actually contribute to order m^a^/A = ma^ (m/A). 
Let us mention that this type of diagrams always contain Coulomb interac- 
tions on external legs and the corresponding integrals are trivial to evaluate. 
All dependence on the infrared cutoff are due to the sensitivity of theory to 
low momenta. Since, in that region, the effective theory is "engineered" to 
reproduce the full theory (here, QED), these singularities are always present 
in the scattering diagrams of both theories so that the renormalized effective 
couplings are always A independent. 

The use of these rules is better explained with the help of a concrete 
example. This is what we do in the next section, where we perform an 
explicit calculation. 

5 Positronium hyperfine splitting 

5.1 Order ma^ 

In this section, we will apply NRQED to a specific calculation, namely the 
hyperfine splitting^ (E (triplet state) - E (singlet state)) of ground state 
positronium. To lowest order, one has the energy levels of the Bohr theory 



^Which we will denote by "hfs" in the following. 
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and there is therefore no hfs. To obtain the first non-zero contribution to the 
sphtting, one has to apply first order perturbation theory to the interactions 
present in the Lagrangian, Eq. ^ . The only interactions that will contribute 
to the order of interest and that are spin-dependent are 

Ci-^o--B, (8) 

2m 

which couples to the transverse photon, and the spin-dependent contact in- 
teraction involving the electron and positron degrees of freedom (or "annihi- 
lation" diagram) 

e2 

4^ X"^cr^ ■ V'^crX- (9) 

m^ 

All the other interactions, because of the presence of additional powers of 

1/m, will not contribute to this order. 

Now we must fix the coefficients Ci and ^3 before computing the hfs. 

The first coefficient, Ci, can be determined by evaluating the spin- flipping 

scattering of an electron from an external field, in the extreme nonrelativistic 

limit. Or one can also simply look up the Pauli Hamiltonian, which the first 

few terms of Eq.0, at tree level, must reproduce. The end result is that 

Ci = 1. The corresponding Feynman rule, in momentum space, is 

± eiP-^Plii^ (10) 

2m 

where ± correspond to the electron and the positron, respectively, and p', p 
are the fermion three-momenta after and before the interaction, respectively. 
The corresponding contribution to the hfs is 



(27r)6 (p2+y)2(q2 + ^2)2 

V''^[^(q - P) X ^]*^ ^^f~2^*^^^ + P) X (^]jX 

^i _ (p-q).(p-q)i ) (11) 

(p-q)2 ^ (p_q)2 

where we have approximated the transverse photon propagator i/{p — qf — 
— z/(p — q)^. Physically, this corresponds to neglecting the recoil since only 
momentum, and no energy, is transferred. The corrections to this approxi- 
mation are of higher order in a. 
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The second term in the square bracket is zero because of the triple prod- 
uct. We can therefore contract the indices i and j. We then have to spin 
average the following expression in the ortho and para states: 

^■f[(p - q) X a]ilj ■ x^[(p - q) X a]x. 

The results are respectively — 2/3(p — q)^ and 2(p — q)^. Taking the first 
value minus the second to get the hfs, we obtain for the integral 

San ffpfq (STrVy/^)^ (p - q)^ _ W 

3^7 "(2^ (p' + 7')2(q2 + 72)Mp-q)' ~ ^~' 
Since the p and q integrals decouple, the calculation is trivial, leading to the 
famous Fermi splitting 0. 

Let us now consider the annihilation interaction. To fix ds, we once more 
turn to the corresponding QED amplitude, namely the tree level annihila- 
tion diagram. In the limit of vanishing external three-momenta, and with 
nonrelativistic normalization for the spinors {u'-jqU = 1 instead of 2E), one 
easily finds that d^ = —ie'^/4. This interaction therefore contributes to the 
hfs 

^^><^ = V 7^^ (p2 + y )2(q2 + ^2y X X ^ XV^ • .P^ax 

= i— X — — xV?/' ■ VVx (13) 

where ■?/' >< ?/) is to be taken as representing the contact interaction sand- 
wiched between two wavefunctions (see the first diagram of fig. [4]). The spin 
average of the spinor expression in Eq. |T^ gives 2 in the ortho state and zero 



in the para state (as it must since, by charge conjugation invariance, para- 
positronium cannot decay to an odd number of photons). The contribution 
from the one annihilation kernel is therefore the well known result: 



ma^ 



+ -^- (14) 



5.2 Order ma^ 



So far, everything has been particularly simple. The evaluation of the a^ 
corrections, however, is a little bit more involved. First, let us notice that, 
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according to the counting rules given above, it is not possible to identify any 
diagram, either in first or second order perturbation theory, that would lead 
to a contribution of order ma^. As an example, consider the annihilation 
interaction in second order perturbation theory {i.e two of these interactions 
sewn together). The wavefunctions would decouple from the inner loop and 
lead to a factor of |\l/(0)p = (ma)^/87r. There would be a factor of (e^/m^)^ 
coming from the two contact interactions, and a factor of m from the nonrel- 
ativistic propagator. In all we obtain an overall factor of a^ times an integral 
in which the only energy scale left is fia. The final result will therefore be 
of order ma^. One can convince oneself that no other diagram can lead to a 
contribution of order ma^. 

Another problem that arises is that, as is well-known from nonrelativistic 
quantum mechanics, divergences arise when we apply second order pertur- 
bation theory to the Lagrangian, Eq.[^. 

These two difficulties are closely related since they have their origins in 
the fact that we have so far entirely neglected relativistic momenta in our 
analysis. This cannot be correct when one goes beyond first order perturba- 
tion theory since loop momenta are allowed to become arbitrarily high. As 
explained above, the correct way to handle this is to renormalize the "bare" 
coefficients of NRQED by matching scattering amplitudes in the effective 
theory with QED diagrams. 

We will now illustrate this procedure in more details with the contribution 
from the annihilation channel, which is represented by Eq. [Q in NRQED. 
Renormalizing the coefficient of the contact interaction to one loop means 
that we impose the sum of all one-loop NRQED diagrams involving this 
interaction to be equal to the corresponding one-loop QED diagrams, as 
shown in fig. [4] . As explained above, the counting rules tell us that there are 
no NRQED bound state diagrams contributing to order ma^. We know that 
the same must hold true for the NRQED scattering diagrams except if there 
are infrared divergent amplitudes. It turns out that there is indeed such an 
amplitude here, given by the annihilation interaction followed by a Coulomb 
exchange. The corresponding integral is equal to 



k2 k2 + A2 

which can easily be seen to be linearly infrared divergent. It will therefore 
contribute to order m'^a^/X and must be included in our calculation. 
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The correction to the lowest order coefficient of the contact interaction 
is then equal to the one-loop QED scattering amplitudes minus the NRQED 
diagram representing the annihilation interaction followed by one Coulomb 
exchange. It is easy to evaluate the QED diagrams because the external 
particles are on-shell and at threshold. The result is 



• 4 -4 

le te m 



27r2m2 47r A 



(15) 



for the vertex interaction, after wavefunction renormalization. Notice the 
linear infrared divergence. It does not appear in the result of the one-loop 
renormalized vertex function T^{p,q) in |]I| (see Eq.[A20]). This is, however, 
due to an error in their calculation as can be seen from the fact that their 
integral Eq.[A17] does contain a linear infrared divergence, divergence that 
cannot be canceled by the vertex renormalization constant which is only 
logarithmically divergent in the infrared. 

For the diagram containing the one-loop vacuum polarization, we find 



• 4 



■" Tn -^ 



97r^m- 



(16) 



The NRQED diagram containing one Coulomb exchange is also easily eval- 
uated and turns out to be 



ie^ m 



(17) 



An X 
so that C(^) is 

.i^(32+y)]/>< (18) 

where >< stands for the (spin averaged) lowest order contact interaction, 
namely —ie^jlrn?. 

The contribution to the hfs is then, up to one loop, equal to 

V'><V^(1 + CW) = 1^(1 + CW) 

ma^ mot' , 2, , . 

= —-—<! + 5' (''" 

which is indeed the correct result 0. 
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5.3 Order ma^ 

We now briefly outline the calculation of the order ma^ correction due to the 
one-photon annihilation diagram. The first thing to do is to renormalize the 
coefficient of the NRQED contact interaction to two loops, as we did at the 
one-loop level in fig. [4]. Keeping only the terms of order a^, we obtain the 
following relation 

>< ^(2) + C^^) X one loop NRQED diagrams 
+ two loops NRQED diagrams 
= QED two loops diagrams (20) 

where all diagrams are scattering amplitudes evaluated at threshold and, 
as mentioned above, we are only considering corrections to the one-photon 
annihilation diagram. One uses once more the counting rules to determine 
which NRQED diagrams will contribute to 0{ma^) in Eq.pO|, taking into 
account the fact that C'-^^ is of order a. 

The rhs (the QED calculation) in Eq. pO[ will be of the form 



2 

A + B ln(A/m) +C- + D— (21) 

A A 

where A, B, C andZ^ are constants of order a^. On the other hand, the 
NRQED diagrams will lead to a result of the form 

2 

A' (A) + B ln(A/7) + C- + D^ (22) 

A A 

where A' {A) contains a finite piece of order ma^ and ultraviolet divergent 
terms. All A-dependent terms are reproduced by the NRQED diagrams so 
that C*^^-* is, in the end, infrared finite (but ultraviolet divergent). The final 
expression for C^'^^ is therefore 

C(2) = {A-A'{A) + B\n{a/2))/>< 

{A-A'{A) + B\n{a/2)). (23) 



ze^ 



Notice that, since the NRQED diagrams contain finite pieces, C-^-* is not 
simply equal to the infrared finite part of Eq. ^1| , A. 
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Once C*^^^ is computed, the correction to the hfs is found by calculating 

iIj><^ (1 + CW + C(2)) 

+ one and two loops NRQED bound state diagrams 

ma^ ma^ , 2, 
= — -— <^+9> 

+ B\n{a/2) + 0{ma^) (24) 

which is now completely ultraviolet finite. It is easy to understand, heuristi- 
cally, why all ultraviolet divergences drop out of Eq.[^. In a nutshell, ultra- 
violet divergent terms arise in NRQED for momenta of order m or greater. 
Since these momenta are much bigger than the typical bound state scale 
7 ~ ma, the UV divergences in the NRQED bound state diagrams can be 
expected to be the same as the ones encountered in the corresponding scatter- 
ing diagrams. Since, in Eq. |^ , each bound state UV divergence is subtracted 



from the corresponding NRQED scattering divergence, all A dependence gets 
canceled^. 

To be more explicit, let us consider in more details one specific NRQED 
diagram. Take, for example, the bound state diagram containing the contact 
interaction iterated twice [i.e. in second order of perturbation theory), as 
depicted in fig. [5]. It is easy to check that the counting rules predict a 
contribution of order ma^ for this diagram. However, the explicit expression 
for this diagram is not well-defined as it diverges linearly: 

d^pd?qd?k 647r7^ e^ 2 



(27r)9 ((k-p)2 + 72)2(q2 + y)2V4^2^ 



1 



—^"^/m — p2/m 



X^a.i) Tr(cricrj) ^VjX 



spinl 



spinO 



ma^ . , ma 



in 



A/m + ^ (25) 



where we have put a cut-off A on the three-momentum flowing through one 
of the contact interaction. This divergence arises from the fact that the 



^Strictly speaking, this is true only for the power-law UV divergences. The heavy scale 
is (implicitly) present in Eq. [24) , where it cancels the m of the 7 factor in the logarithm. 
See our previous footnote on the running of the low energy coefficients. For a discussion 
of the issue of power-law divergences in a broader context, see jj. 
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Lagrangian Eq. [0] is defined for momenta mucli smaller than the electron 
mass whereas the momentum in a loop is allowed to become relativistic. 
However, for p > m, one expects the intermediate state loop to shrink to 
a point (in coordinate space), as far as the low energy theory is concerned. 
This means that the divergent contributions to the NRQED bound state 
diagrams can be canceled by a renormalization of the tree level coefficients. 
In fact, the necessary counterterms are already present in the coefficients 
C^^^ and C^'^\ To see this, let us go back to the definition of C^'^\ fig-[4]- 
When we calculated the ma^ correction, the only one-loop NRQED diagram 
we kept in this equation was the contact interaction followed by a Coulomb 
exchange. Now, however, in order to be consistent in our calculation of order 
ma^, we must go back to that equation and include all necessary NRQED 
one-loop diagrams. One of those diagrams turns out to be the scattering 
diagram containing the contact interaction iterated twice {i.e fig. [5] with free 
spinors on the external legs instead of wavefunctions), which has the simple 
Feynman rule 

(26) 



/ 



(27r)'^ k^ Am? m? m 

This leads to a new contribution to C*-^-' and, consequently, to a correction 
to bound state energy equal to (see Eq.p4|) 



|xl>(0)|2 X -2—- = . 27 

m^ m An m 



Combining this with Eq.|25] yields a finite correction to the hyperfine split- 
ting equal to ma^/S. 

Before wrapping up this section, we would like to mention that there is 
yet another new ingredient appearing at order ma^. It has to do with the 
fact that, as we saw earlier, the Coulomb interaction does not increase the 
order (in powers of a) of the bound state diagrams. Although, in first order 
perturbation theory this is taken care of by the Schrodinger wavefunctions, 
in second order one must also sum up the Coulomb interactions in the inter- 
mediate state. In practice, this is done by using a closed form expression for 
the full Coulomb propagator which has been derived in both coordinate and 
momentum space H]. 



6 Discussion 

We have seen that there are two classes of contributions to bound state 
properties, one being characterized by the typical bound state momentum 
p ?a ma and one characterized by relativistic momenta p ^ m. Such a 
situation is ideal for the use of an effective field theory (in the decoupling 
scenario). NRQED is such a theory, with heavy scale A equal to the electron 
mass m. The use of NRQED greatly simplifies bound state calculations since 
it involves only nonrelativistic interactions in the bound state diagrams and 
incorporates the relativistic physics via the matching of conventional scat- 
tering amplitudes. The relative simplicity of the nonrelativistic interactions 
makes the bound state diagrams easy to evaluate whereas the QED inter- 
actions appear only in scattering amplitudes, which can be evaluated using 
conventional methods. This is to be contrasted with traditional approaches, 
the Bethe-Salpeter equation for example, which involve the fully relativistic 
theory in the bound state diagrams, leading to the difficulties discussed in 
Section 2. 

We are finally able to understand in what way the order ma'^ and ma^ 
corrections have a special status: the ma^ corrections come exclusively from 
nonrelativistic momenta {i.e. from NRQED diagrams) whereas the ma^ 
corrections come practically exclusively from relativistic momenta [i.e. from 
QED scattering diagrams). We say "practically" because, as we saw in the 
example derived above, there is a linear infrared divergence in the QED 
diagrams, showing that they are in fact sensitive to very low momenta. In 
fact this divergence is a pure threshold singularity, which is due to tha fact 
that the diagrams are evaluated with external three-momenta set to zero. 
However, being a low-energy phenomenon, one expects this divergence to be 
reproduced by NRQED and we saw that this is exactly what happens, in 
the form of the Coulomb exchange diagram. But aside from this singularity, 
which gets canceled when matching QED and NRQED, the ma^ contribution 
originates entirely from QED scattering diagramsQ Q 

^"^ Although we have explicitly shown this only for the annihilation diagram, this also 
holds for the p ■ A interaction; therefore all 0{mQ^) contributions to the hfs are purely 
relativistic. 

^^It is in that sense that we call these contributions "purely" relativistic. Maybe a 
more appropriate terminology would be "bound state independent" . Notice that what we 
call "relativistic contributions" are sometimes referred to as "radiative corrections" in the 
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Let us now compare to the technique of Ref. |^. For the 0{ma'^) cor- 
rection, they consider the QED one-photon exchange in the scattering and 
annihilation channels, expanded them to lowest order in p/m, and obtained 
the expressions Eq.|jlT|] and Eq.[0]. At this order, this is exactly equivalent 



to using NRQED. In the language of NRQED, these interactions are found 
in the effective Lagrangian Eq. [Q, which is obtained by taking the nonrel- 
ativistic limit of the QED Lagrangian. It is in that sense that we call the 
0{ma'^) "purely" nonrelativistic; they are found directly from the effective 
Lagrangian. 

Now consider the 0{ma^) corrections. Once again, we restrict ourselves 
to the one-photon annihilation corrections for the sake of simplicity. In ref. [||], 
these corrections are obtained by sandwiching QED one-loop diagrams be- 
tween wavefunctions, and taking the limit p — > in the QED spinors (thresh- 
old limit). In that limit, their calculation reduces to multiplying the one- loop 
threshold QED diagrams by the square of the wavefunctions evaluated at 
the origin. This is almost exactly what we find using NRQED. Indeed the 
{ma^) term comes uniquely from the renormalization of the NRQED "bare" 
parameters, as we demonstrated explicitly in this paper. However, there is an 
important distinction between the two calculations. In |]l|, the photon mass 
dependence was discarded, with no further justification. But in NRQED, the 
A dependence gets canceled, order by order in the matching procedure. Also, 
it is now easy to understand why the 0{ma^) result for the ground state 
can be trivially generalized to ( / = 0) states of arbitrary principal quantum 
number. Indeed, since the only effect of the bound state is to provide an 
overall normalization equal to the square of the wavefunction at the origin, 
the 0{ma^) will carry the same dependence on n as |\E'(0)p, namely it is 
proportional to 1/n^. 

We are now in position to understand how the two approaches would dif- 
fer at 0{ma^), not only in the handling of the infrared divergences but also in 



the prediction of the finite term. We saw in Eq.|]2^ that, in contradistinction 
to the 0{ma^) and 0{ma^) corrections, the ma^ term contains finite con- 
tributions from both relativistic and nonrelativistic momenta. Notice that 
this is not simply a consequence of the presence of loops since the 0{ma^) 
term already included integration over internal momenta. It has for origin 
a subtle interplay between relativistic and nonrelativistic flow of momenta 

literature. 
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through the Feynman diagrams, interplay which is brought in evidence by the 
use of NRQED and its counting rules. Using the rules of ^ would yield the 
relativistic finite contribution only (in addition to infrared divergent terms 
which would have to be discarded without further justification). 

The fact that nonrelativistic momenta play an important role at this or- 
der also has for consequence that it is no longer possible to consider Feynman 
diagrams containing only a fixed number of loops, as done in [^ . Indeed, as 
already mentioned in Section 5.3, at this order one must consider bound state 
diagrams with an infinite number of Coulomb interactions in the intermedi- 
ate states. For example, were we to complete our calculation of the 0{ma^) 
contribution due to the annihilation interaction in second order of perturba- 
tion theory, we would have to sum up the same diagram with any number 
of Coulomb vertices between the two contact interactions. This complica- 
tion is directly at the origin of the issue of reducible vs irreducible Feynman 
diagrams in the Bethe-Salpeter approach. In NRQED, with the help of the 
Coulomb gauge, this problem is reduced to its most simple form since it is 
confined to the Coulomb interaction and the infinite sum can be performed 
using the expressions of p. 

Addendum 

While writing this paper we became aware of 0, in which one contribution of 
0{ma^) is calculated. This contribution corresponds to the one-loop vacuum 
polarization correction to the exchange of two photons. 

These diagrams, however, lead to purely relativistic corrections or, in 
other words, they contribute only to the coefficient A in Eq. ||2^. This is 
because vacuum polarization is a highly virtual effect (notice that there is 
no vacuum polarization per se in NRQED; it only enters through the renor- 
malization of the low energy theory's coefficients). This can be seen directly 
from the fact that the corresponding correction to the bound state energy is 
proportional to the square of the wavefunction at the origin (see Eqs[14-15] 
in ; notice that K is independent of the bound state momenta p and q) . In 
that respect, the corrections due to these diagrams are similar to the 0{ma^) 
terms; they are insensitive to bound state physics and can be treated using 
the formalism of |1|] (especially since they are infrared finite). An example of 
a diagram contributing to 0{ma^) that would be sensitive to low momenta 
would be, for example, the three photon exchange diagrams. 
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